LA-UR 09-06331 
INT-PUB-10-054 



The Electrosphere of Macroscopic "Quark Nuclei": 
A Source for Diffuse MeV Emissions from Dark Matter. 

Michael McNeil Forbes* 

Institute for Nuclear Theory, University of Washington, Box 351550, Seattle, Washington, 98195-1550, USA 

Kyle Lawson and Ariel R. Zhitnitsky 

Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, V6T 1Z1, Canada 

(Dated: October 5, 2010) 

Using a Thomas-Fermi model, we calculate the structure of the electrosphere of the quark antimat- 
ter nuggets postulated to comprise much of the dark matter. This provides a single self-consistent 
density profile from ultrarelativistic densities to the nonrelativistic Boltzmann regime that use to 
present microscopically justified calculations of several properties of the nuggets, including their 
net charge, and the ratio of MeV to 511 keV emissions from electron annihilation. We find that the 
calculated parameters agree with previous phenomenological estimates based on the observational 
supposition that the nuggets are a source of several unexplained diffuse emissions from the Galaxy. 
As no phenomenological parameters are required to describe these observations, the calculation pro- 
vides another nontrivial verification of the dark-matter proposal. The structure of the electrosphere 
is quite general and will also be valid at the surface of strange-quark stars, should they exist. 
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I. INTRODUCTION 

In this paper we explore some details of a testable and 
well-constrained model for dark matter [1-5] in the 
form of quark matter as antimatter nuggets. In par- 
ticular, we focus on physics of the "electrosphere" sur- 
rounding these nuggets: It is from here that observable 
emissions emanate, allowing for the direct detection of 
these dark-antimatter nuggets. 

We first provide a brief review of our proposal in 
Sec. II, then describe the structure of the electrosphere 
of the nuggets using a Thomas-Fermi model in Sec. IV. 
This allows us to calculate the charge of the nuggets, 
and to discuss how they maintain charge equilibrium 
with the environment. We then apply these results to 
the calculation of emissions from electron annihilation 
in Sec. V, computing some of the phenomenological pa- 
rameters introduced in [6] required to explain current 
observations. The values computed in the present pa- 
per are consistent with these phenomenologically moti- 
vated values, providing further validation of our model 
for dark matter. The present results concerning the den- 
sity profile of the electrosphere may also play an im- 



2 



portant role in the study of the surface of quark stars, 
should they exist. 

II. DARK MATTER AS DENSE QUARK NUGGETS 

Two of the outstanding cosmological mysteries - the 
natures of dark matter and baryogenesis - might be ex- 
plained by the idea that dark matter consists of com- 
pact composite objects (ccos) [1-5] similar to Witten's 
strangelets [7]. The basic idea is that these ccos - 
nuggets of dense matter and antimatter - form at the 
same qcd phase transition as conventional baryons 
(neutrons and protons), providing a natural explana- 
tion for the similar scales Q DM ~ 5Qg. Baryogenesis 
proceeds through a charge separation mechanism: both 
matter and antimatter nuggets form, but the natural 
cp violation of the so-called 6 term in qcd 1 - which 
was of order unity 6 ~ 1 during the qcd phase transi- 
tion - drives the formation of more antimatter nuggets 
than matter nuggets, resulting in the leftover baryonic 
matter that forms visible matter today (see [2] for de- 
tails). Note, it is crucial for our mechanism that cp vio- 
lation be able to drive charge separation: though not yet 
proven, this idea may already have found experimen- 
tal support through the Relativistic Heavy Ion Collider 
(rhic) at Brookhaven [13], where charge separation ef- 
fects seem to have been observed [14, 15] 

The mechanism requires no fundamental baryon 
asymmetry to explain the observed matter /antimatter 
asymmetry. Together with the observed relation Q DM « 
5Qg (see [16] for a review) we have 

^universe = = £> nU gg e t + ^visible — ^antinugget (^ a ) 
^dark-matter = ^nugget + S an ti nu gg e t ~ ^visible (^ D ) 

where B universe is the overall asymmetry - the total 
number of baryons 2 minus the number of antibaryons 
in the Universe - and £>dark-matter is the total number 
of baryons plus antibaryons hidden in the dark-matter 



1 If 6 is nonzero, one must confront the so-called strong cp problem 
whereby some mechanism must be found to make the effective 9 
parameter extremely small today in accordance with measurements. 
This problem remains one of the most outstanding puzzles of the 
Standard Model, and one of the most natural resolutions is to in- 
troduce an axion field. (See the original papers [8-10], and recent 
reviews [11].) Axion domain walls associated with this field (or 
ultimately, whatever mechanism resolves the strong cp problem) 
play an important role in forming these nuggets, and may play in 
important role in their ultimate stability. See [1, 2, 12] for details. 

2 Note that we use the term "baryon" to refer in general to anything 
carrying Ug(l) baryonic charge. This includes conventional colour 
singlet hadrons such as protons and neutrons, but also includes 
their constituents - i.e. the quarks - in other phases such as strange 
quark matter. 



nuggets. The dark matter comprises a baryon charge 
of £> nu gget contained in matter nuggets, and an an- 
tibaryonic charge of S ant i nu gg et contained in antimatter 
nuggets. The remaining unconfined charge of B v isible is 
the residual 'Visible" baryon excess that forms the reg- 
ular matter in our Universe today. Solving Eq. (1) gives 
the approximate ratios Bantinugget^nugget ^visible - 3:2:1. 

Unlike conventional dark-matter candidates, dark- 
matter /antimatter nuggets will be strongly interacting, 
but macroscopically large, objects. They do not contra- 
dict any of the many known observational constraints 
on dark matter or antimatter [3] for three reasons: 

1. They carry a huge (anti)baryon charge \B\ « 10 20 - 10 30 , 
so they have an extremely tiny number density. 3 This 
explains why they have not been directly observed on 
earth. The local number density of dark-matter parti- 
cles with these masses is small enough that interactions 
with detectors are exceedingly rare and fall within all 
known detector and seismic constraints [3]. (See also 
[17, 18] and references therein.) 

2. The nugget cores are a few times nuclear density p ~ 
10 GeV/fm 3 , and thus have a size R ~ 10 -7 -10 -3 cm. 
Their interaction cross section is thus small a/M ~ 
4:7tR 2 /M = 10~ 13 -10~ 9 cm 2 /g: well below the typical 
astrophysical and cosmological limits, which are on the 
order of a/M < 1 cm 2 /g. Dark-matter-dark-matter in- 
teractions between these nuggets are thus negligible. 

3. They have a large binding energy such that the bary- 
onic matter in the nuggets is not available to participate 
in big bang nucleosynthesis (bbn) at T « 1 MeV. In par- 
ticular, we suspect that the core of the nuggets forms a 
superfluid with a gap of the order A ^ 100 MeV, and 
critical temperature T c ~ A/ y/2 « 60 MeV, as this scale 
provides a natural explanation for the observed photon 
to baryon ratio ftg/ft 7 ~ 10 -10 [2], which requires a 
formation temperature of Tf orm = 41 MeV [19]. 4 

Thus, on large scales, the nuggets are sufficiently 
dilute that they behave as standard collisionless cold 
dark matter (ccdm). When the number densities of 
both dark and visible matter become sufficiently high, 
however, dark-antimatter-visible-matter collisions may 
release significant radiation and energy. In particular, 
antimatter nuggets provide a site at which interstellar 
baryonic matter - mostly hydrogen - can annihilate, 
producing emissions that should be observable from 
the core of our Galaxy of calculable spectra and energy. 



3 If the average nugget size ends up in the lower range, then the 
Pierre Auger observatory may provide an ideal venue for searching 
for these dark-matter candidates. 

4 At temperatures below the gap, incident baryons with energies be- 
low the gap would Andreev reflect rather than become incorpo- 
rated into the nugget. 
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These emissions are not only consistent with current 
observations, but naturally explain several mysterious 
diffuse emissions observed from the core of our Galaxy, 
with frequencies ranging over 12 orders of magnitude. 

Although somewhat unconventional, this idea natu- 
rally explains several coincidences, is consistent with 
all known cosmological constraints, and makes testable 
predictions. Furthermore, this idea is almost entirely 
rooted in conventional and well-established physics. In 
particular, there are no "free parameters" that can be 
- or need to be - "tuned" to explain observations: In 
principle, everything is calculable from well-established 
properties of qcd and qed. In practice, fully calculat- 
ing the properties of these nuggets requires solving the 
fermion many-body problem at strong coupling, so we 
have generally resorted to "fitting" a handful of phe- 
nomenological parameters from observations. In this 
paper we examine the qed physics of the electrosphere, 
providing a microscopic basis for some of these parame- 
ters. Once these parameters are determined, the model 
makes unambiguous predictions about other processes 
ranging over more than 10 orders of magnitude in scale. 

The basic picture involves the antimatter nuggets - 
compact cores of nuclear or strange-quark matter (see 
Sec. Ill) surrounded by a positron cloud with a profile 
as calculated in Sec. IV. Incident matter will annihilate 
on these nuggets producing radiation at a rate propor- 
tional to the annihilation rate, thus scaling as the prod- 
uct pv(r)p BM (r) of the local visible and dark-matter 
densities. This will be greatest in the core of the Galaxy. 
To date, we have considered five independent observa- 
tions of diffuse radiation from the core of our Galaxy: 

1. Spi/ integral (the spi instrument on the integral satel- 
lite) observes 511 keV photons from positronium decay 
that is difficult to explain with conventional astrophys- 
ical positron sources [20-22]. Dark-antimatter nuggets 
would provide an unlimited source of positrons as sug- 
gested in [23, 24]. 

2. Comptel/crgo (the comptel instrument on the Comp- 
ton Gamma Ray Observatory (crgo) satellite) detects 
a puzzling excess of 1-20 MeV 7-ray radiation. It was 
shown in [6] that the direct e + e~ annihilation spectrum 
could nicely explain this deficit, but the annihilation 
rates were crudely estimated in terms of some phe- 
nomenological parameters. In this paper we provide 
a microscopic calculation of these parameters (Eqs. (19) 
and (21)), thereby validating this prediction. 

3. Chandra (the Chandra x-ray observatory) observes a dif- 
fuse keV x-ray emission that greatly exceeds the en- 
ergy from identified sources [25]. Visible-matter /dark- 
antimatter annihilation would provide this energy. It 
was shown in [4] that the intensity of this emission is 
consistent with the 511 keV emission if the rate of pro- 
ton annihilation is slightly suppressed relative to the 
rate of electron annihilation. In Sec. IV C we describe 



the microscopic nature of this suppression. 

4. Egret /crgo (the Energetic Gamma Ray Experiment 
Telescope aboard the crgo satellite) detects MeV to GeV 
gamma rays, constraining antimatter annihilation rates. 
It was shown in [4] that these constraints are consistent 
with the rates inferred from the other emissions. 

5. Wmap (the Wilkinson Microwave Anisotropy Probe) 
has detected an excess of GHz microwave radiation 
- dubbed the "wmap haze" - from the inner 20° 
core of our Galaxy [26-29]. Annihilation energy not 
immediately released by the above mechanisms will 
thermalize, and subsequently be released as thermal 
bremsstrahlung emission at the eV scale. In [5] it was 
shown that the predicted emission from the antimatter 
nuggets is consistent with, and could completely ex- 
plain, the observed wmap haze. 

These emissions arise from the following mechanism: 
Neutral hydrogen from the interstellar medium (ism) 
will easily penetrate into the electrosphere, providing a 
source of electrons and protons. 

The first and simplest process is the annihilation of 
the electrons through positronium formation, produc- 
ing 511 keV photons as discussed in [23, 24]. Note that 
this mechanism predicts that, within the environment 
of the electrosphere, virtually all of the low-energy 
emission should be characterized by a positronium de- 
cay spectrum, including 25% as a sharp 511 keV line 
from the decay of the singlet state and the remaining 
75% as the broad three photon continuum resulting 
from the triplet state. As emphasized in [30], simply 
postulating a dark-matter source of positrons does not 
suffice to explain the observed spectrum characterized 
by 94 ± 4% positronium annihilation: The positrons 
must annihilate in the appropriate cool environment as 
provided by the electrosphere of the nuggets. 

Our proposal can thus easily explain the observed 
511 keV radiation. If we assume that this process is the 
dominant source, then we can use this to normalize the 
intensities of the other emissions. A remarkable feature 
of this proposal is that it then predicts the correct inten- 
sity for all of the other observations, even though they 
span many orders of magnitude in frequency. 

The second process is direct annihilation of the 
electrons on the positrons in the electrosphere. As 
discussed in appendix B, the electrons are strongly 
screened by the positron background, and some frac- 
tion can penetrate deep within the electrosphere. There 
they can directly annihilate with high-momentum 
positrons in the Fermi sea producing radiation up to 
10 MeV or so. This process was originally discussed 
in [6] where the ratio of direct MeV annihilation to 
511 keV annihilation was characterized by several phe- 
nomenological parameters chosen to fit the observa- 
tions. In Sec. V we put this prediction on solid ground 
and show that these parameter fits agree with the micro- 
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scopic calculation based on the electrosphere structure, 
which depends only on qed. 

The other radiation originates from the energy de- 
posited by annihilation of the incident protons in the 
core of the nuggets. In [4] we argue that the protons 
will annihilate just inside the surface of the core, releas- 
ing some 2 GeV of energy Occasionally this process 
will release GeV photons - the rate of which is consis- 
tent with the egret/ crgo constraints - but most of the 
energy will be transferred to strongly interacting com- 
ponents, and ultimately about half will scatter down 
into ~ 5 MeV positrons that stream out of the core. (The 
5 MeV scale comes from the effective mass of the photon 
in the medium which mediates the energy exchange). 
These positrons are accelerated by the strong electric 
fields, and emit field-induced bremsstrahlung emission 
in the 10 keV band - a scale set by balancing the rate 
of emission with the local plasma frequency (photons 
can only be emitted once the plasma frequency is low 
enough). The spectrum is also calculable: It is very flat 
and similar to a thermal spectrum without the sharp 
falloff above the temperature scale. This is consistent 
with the Chandra observations which cannot resolve the 
thermal falloff, but future analysis might be able to dis- 
tinguish between the two. 

We shall shed some light on the interaction between 
the protons and the core in Sec. IV C, but cannot yet 
perform the required many-body analysis to place all of 
this on a strong footing as this would require a practical 
solution to outstanding problems of high-density qcd. 

The remaining energy will thermalize within the 
nuggets, until an equilibrium temperature of about 
T ~ 1 eV is reached [5]. Direct thermal emission at 
this scale would be virtually impossible to see against 
the backgrounds, but the spectrum - calculable entirely 
from qed - extends to very low frequencies, and the 
intensity of the emission in the microwave band is just 
enough to explain the wmap haze [5]. 



III. CORE STRUCTURE 

A full accounting of this nugget proposal requires a 
proper description of the high-density phase found in 
the core of the nuggets. Unfortunately, a quantitative 
understanding of this phase requires a practical solu- 
tion to the notoriously difficult problem of high-density 
qcd. The density of the core will be within several or- 
ders of nuclear matter density ~ 1 — 100 GeV/ fm 3 . This 
is not high enough for the asymptotic freedom of qcd 
may be used to solve the problem perturbatively, so one 
must resort to nonperturbative techniques such as the 
lattice formulation of qcd. Unfortunately, at finite den- 
sity, the presence of the infamous sign problem renders 
this approach exponentially expensive and it remains a 



famously intractable problem. As such, we cannot ex- 
actly quantify the nature of the core and must constrain 
its properties from other observations; we list the im- 
portant properties in this section. 

Fortunately, the observable emissions discussed in 
this paper result primarily from the calculable physical 
processes in the electrosphere of the nuggets, and are 
thus largely insensitive to the exact nature of the core. 
Nevertheless, the core structure must be addressed, and 
it is possible that future developments concerning the 
properties of high-density qcd could rule out the feasi- 
bility of our nugget proposal. 

The first problem concerns the stability of the core. 
All evidence suggests that, in the absence of an external 
potential (such as the gravitational well of a neutron 
star), nuclear matter will fragment into small nuclei 
with a baryon number no larger than a few hundred. 
This suggests one of two possibilities: 

• The first possibility suggested by Witten [7] is that a 
phase of strange-quark matter [31] becomes stable at 
high density. In this case, the nuggets are simply 
strangelets and antistrangelets and the novel feature 
here is that the domain walls associated with strong 
CP violation provide the required mechanism to con- 
dense enough matter to catalyze the formation of the 
strangelets before they evaporate. (For a brief review, 
see [32] and references therein.) 

Although strangelets have not yet been observed, the 
possibility of stable strange-quark mater has not yet 
been ruled out (see for example [33]). This must be care- 
fully reconciled with future astrophysical observations 
and constraints as it is conceivable that this possibility 
might be ruled out in the future. If the nuggets are a 
form of strange-quark matter, one must also consider 
the possibility of mixed phases as suggested in [34] (see 
for example [35, 36] and references therein). This would 
most likely have to be ruled out energetically to prevent 
the nuggets from fragmenting. 

• The second possibility is that the domain walls responsi- 
ble for forming the nuggets at the qcd phase transition 
become an integral part, providing a surface tension 
that holds the nuggets together, even in the absence of 
absolutely stable strange-quark matter. The stability of 
this possibility has been discussed in detail in [1] and 
we shall not repeat these arguments here. In this case, 
the core may be something more akin to dense nuclear 
matter as might be found in the core of a neutron star. 

As discussed above, in order to explain the observa- 
tions, the following core properties are crucial to our 
proposal. This provides some insight into the required 
nature of the core: 

1. The nuggets must be stable. Stable strange-quark mat- 
ter would offer a nice explanation. Otherwise, the struc- 
ture of the core - especially the surface - must be con- 
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sidered in more detail. This is a complicated problem 
that we do not presently know how to solve. 

2. As discussed above, the formation of the objects must 
stop at T = 41 MeV to explain the observed photon 
to baryon ratio. This could be naturally explained by 
the order 100 MeV pairing gaps expected in colour- 
superconducting strange-quark matter. If the proposal 
turns out to be correct, then, this would provide a 
precise measurement of the pairing properties of high- 
density qcd. (The exact relationship between the pair- 
ing gap and the formation temperature will be quite 
nontrivial and require a detailed model of the forma- 
tion dynamics.) This favours a model of the core with 
strong pairing correlations. 

3. In order to explain the Chandra data, our picture of the 
emission mechanism requires that the proton annihila- 
tions occur somewhat within the core, not immediately 
on the surface where there would be copious pion pro- 
duction. The simplest explanation for this would be the 
presence of strong correlations in the core, delaying the 
annihilation until the protons penetrate a few hundred 
fm or so. Again, if the core is a colour superconductor, 
then the pairing correlations could explain this, but a 
detailed calculation is needed to make sure. 

An argument supporting the presence of such correla- 
tions in a dense environment is the observation that the 
annihilation cross section of an antiproton on a heavy 
nuclei is many times smaller than the vacuum pp and 
nn annihilation rates suggest. Indeed, it is expected that 
antiproton /nuclei bound states may persist with a life- 
time much longer than the vacuum annihilation rates 
predicts (see [37, 38]). This annihilation suppression 
should become even stronger with increasing density. 

Thus, the hypothesis of strange-quark matter simplifies 
the picture quite a bit, but we do not yet see a way to di- 
rectly link this hypothesis with the existence or stability 
of the nuggets. 



IV. ELECTROSPHERE STRUCTURE 

Here we discuss the density profile of the electron cloud 
surrounding extremely heavy macroscopic "nuclei" We 
present a type of Thomas-Fermi analysis including the 
full relativistic electron equation of state required to 
model the relativistic regime close to the nugget core. 
We consider here the limit when the temperature is 
much smaller than the mass T <C 511 MeV. The solu- 
tion for higher temperatures follows from similar tech- 
niques with fewer complications. 

The observable properties of the antimatter nuggets 
discussed so far [1-6] depend on the existence of a non- 
relativistic "Boltzmann" regime with a density depen- 
dence n(r) ~ (r — tb) -2 (see appendix 1 of Ref. [5] 
for details). This region plays an important role in ex- 



plaining the wmap haze [5] as well as in the analysis of 
the diffuse 511 keV emissions [2]. The techniques pre- 
viously used, however, were not sufficient to connect 
this nonrelativistic "Boltzmann" regime to the relativis- 
tic regime through a self -consistent solution determined 
by parameters T, ]i, m e . 

The main point of this section is to put the exis- 
tence of a sizable Boltzmann regime on a strong foot- 
ing, and to calculate some of the previously estimated 
phenomenological parameters that depend sensitively 
on density profile. These may now be explicitly com- 
puted from qed using a justified Thomas-Fermi approx- 
imation to reliably account for the many-body physics. 
We show that a sizable Boltzmann regime exists for all 
but the smallest nuggets, which are ruled out by lack 
of terrestrial detection observation. We also address 
the question of how the nuggets achieve charge equilib- 
rium (Sec. IV C), discussing briefly the charge-exchange 
mechanism, and determining the overall charge of the 
nuggets (see Table I). 



A. Thomas-Fermi Model 

To model the density profile of the electrosphere, we 
use a Thomas-Fermi model for a Coulomb gas of 
positrons. This is derivable from a density functional 
theory (see Appendix A) after neglecting the exchange 
contribution, which is suppressed by the weak coupling 
oc. The electrostatic potential <p(r) must satisfy the Pois- 
son equation 



\/ 2 (p{r) = -4/ren(r). 



(2) 



where en(r) is the charge density. Outside of the nugget 
core, we express everything in terms of the local effec- 
tive chemical potential 



fi(r) -e<p(r), 



(3) 



and express the local charge density through the func- 
tion q(r) = en[ji(r)] where n\]\\ contains all of the in- 
formation about the equation of state. 

As with the Thomas-Fermi model of an atom, the self- 
consistent solution will be determined by the charge 
density of the core ( // nucleus ,/ ). We shall simply imple- 
ment this as a boundary condition at the nugget core 
boundary at radius r = R, and thus only consider the 
region r > R. The resulting solution may thus be ex- 
pressed in terms of the equations 



V 2 }i(r) = 4tnocn[}i(r)], 



\p\\ 2 



n\p]=2j- 



d 3 p 

(2nf Ll + e ( e p-^) /T l + e ( e p+^ /r . 



(4a) 
, (4b) 



with the appropriate boundary conditions at r = R and 
r = oo. In (4b) we have explicitly included both particle 
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and antiparticle contributions, as well as the spin de- 
generacy factor, and have used modified Planck units 
where h = c = 4zreo = 1 so that e 2 = oc and energy, 
momentum, inverse distance, and inverse time are ex- 
pressed in eV. 

We assume spherical symmetry so that we may write 

v2 FW = ^/Tr^ = f{r)+2Ar)/r. (5) 

Close to the surface of the nuggets, the radius is 
sufficiently large compared with the relevant length 
scales that the curvature term 2ji f (r)/r may be ne- 
glected. We shall call the resulting approximation "one- 
dimensional" The resulting profile does not depend 
on the size of the nuggets and remains valid until the 
electrosphere extends to a distance comparable to the 
radius of the nugget, at which point the full three- 
dimensional form will cut off the density When ap- 
plied to the context of strange stars [39-42], the one- 
dimensional approximation will be completely suffi- 
cient. See also [43] where a Thomas-Fermi calculation 
is used to determine the charge distribution inside a 
strange-quark nugget. 

To help reason about the three-dimensional equation, 
we note that (4a) may be expressed as 

]i"{x) = x~ A 4jtocn[}i(x)] (6) 

where x = 1/r. A nice property of this transformation 
is that ]i(x) must be a convex function if the charge 
density has the same sign everywhere. 

1. Boundary Conditions 

The physical boundary condition at the origin follows 
from smoothness of the potential. By combining this 
with a model of the charge distribution in the core and 
an appropriate long-distance boundary-condition, one 
could in principle model the entire distribution of elec- 
trons throughout the nugget. 

In the nugget core, however, beta-equilibrium essen- 
tially establishes a chemical potential on the order of 10 
- 100 MeV [39] that depends slightly on the exact equa- 
tion of state for the quark-matter phase (which is not 
known). Thus, we may simply take the boundary con- 
dition as }i(R) = }Ir ~ 25 MeV. Our results are not very 
sensitive to the exact value, though if less that 20 MeV, 
then this acts as a cutoff for the direct e + e~ emission 
discussed in Sec. V A 2. 

The formal difficulty in this problem is properly for- 
mulating the long-distance boundary conditions. At 
T = 0, the large-distance boundary condition for non- 
relativistic systems is clear: n(r —> oo) =0. From (6) we 
see that jt(x) is linear with slope y!{x) = — eQ where 



Q is the overall charge of the system (ions with a defi- 
ciency of electrons /positrons are permitted in the the- 
ory; see for example [44]): 

eQ(r) = J df 4:7tocn(f)f 2 = J df (r 2 }i'(f)^ =r 2 ji'(r). 

At finite temperature, however, this type of boundary- 
condition is not appropriate. Instead, one must con- 
sider how equilibrium is established. 

In a true vacuum, a finite temperature nugget will 
"radiate" the loosely bound outer electrosphere until 
the electrostatic potential is comparable to the temper- 
ature eQ/r* rsj T. At this radius, radiation becomes 
exponentially suppressed. Suppose that the density at 
this radius is n(r*). We can estimate an upper bound 
for the evaporation rate as Arcr 2 vn where v ~ y/T/m. 5 

This needs to be compensated by rate of charge de- 
position from the surrounding plasma which can be es- 
timated as 47r# 2 7J N+/ISM ft N+/ISM where ^ N + /ISM ~ 10~ 3 c 
is the typical relative speed of the nuggets and charged 
components in the Inter-Stellar Medium (ism) and 4nR 2 
is approximately the cross section for annihilation on 
the core (see Sec. IV C). The density n N + ISM of charged 
components in the core of the Galaxy is typically 10 _1 
to 10 -2 of the total density n 1SM ~ 1 cm -3 . Thus, once 
the density falls below 

« < n md ~ ~f2 Y Y^ism^ismIO (7) 

charge equilibrium can easily be established. This al- 
lows us to formulate the long-distance boundary condi- 
tions by picking a charge Q and outer radius r* such 
that eQ/r* ~ T, and n(r*) ~ n ra( j, establishing both 
the typical charge of the nuggets as well as the outer 
boundary condition for the differential equation (4). 6 

Note that the Thomas Fermi approximation is not 
trust-worthy in these regimes of extremely low density, 
but the boundary condition suffices to provide an es- 
timate of the charge of the nuggets: were they to be 
less highly charged, then the density at r* would be 
sufficiently high that evaporation would increase the 
charge; were they more highly charged, then the evap- 
oration rate would be exponentially suppressed, allow- 
ing charge to accumulate. 



5 The Boltzmann averaged velocity is lower (v) — \J~T~I (Inrri), and 
cutting off the integral properly will lower this even more, so this 
gives a conservative upper bound on the evaporation rate. 

6 Strictly speaking, at finite T, the equations do not have a formal 
solution with a precise total charge because there is always some 
density for ji > 0. Practically, once m — ji <C T, the density becomes 
exponentially small and the charge is effectively fixed. 
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Thus, for the antimatter nuggets, the density profile is 
determined by system (4) and the boundary conditions 



eQ/r, 



p(R) = jt R ~25MeV, (8a) 
T-leV, n(u) = n rad (r*). (8b) 



where n[ji] oc e^ T : 

«b(z) 



1 



(10) 



2tta: (z + z b ) 2 ' 
Since we have the full density profiles, we fit Zb here 



B. Profiles 

In principle, one should also allow the temperature to 
vary However, tThe rate of radiation in the Boltzmann 
regime is suppressed by almost six orders of magnitude 
with respect to the black-body rate; meanwhile, the den- 
sity of the plasma ~ (mT) 3 ^ 2 is quite large [5]. Thus, 
the abundance of excitations ensures that the thermal 
conductivity is high enough to maintain an essentially 
constant temperature throughout the electrosphere. 

One can now numerically solve the system (4) 
and (8). We shall consider six cases: total (anti)baryon 
charge B G {10 20 , 10 24 , 10 33 } and quark-matter densi- 
ties ft CO re £ {1, 100}fm -3 . The resulting profiles are 
plotted in Fig. I. 7 As a reference, we also plot the "one- 
dimensional" approximation obtained numerically by 
neglecting the curvature term in (5) as well as two 
analytic approximations: One for the ultrarelativistic 
regime [39, 40, 46] 8 where n[ji] « fi 3 /3n 2 : 
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and one for the nonrelativistic Boltzmann regime [5] 
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TABLE I. Antimatter nugget properties: the core radius R and 
total electric charge Q for nuggets with two different assumed 
core densities n cor e and three different antibaryon charges B. 
These correspond to the profiles shown in Fig. 1. 



Note that the ultrarelativistic approximation employed in [45] dis- 
cussing nonrelativistic physics not only overestimates the density 
by 3 orders of magnitude in the Boltzmann regime, but also has 
a different z-dependence [(9) verses (10)]: The ultrarelativistic ap- 
proximation is not valid for ji < m where most of the relevant 
physical processes take place. 
For T 0, see appendix A 2 a. 



FIG. 1. Density profile of antimatter nuggets (bottom plot is 
a zoom). The thick solid line is the one-dimensional approxi- 
mation neglecting the curvature of the nugget. Six profiles are 
shown descending from this in pairs of thin solid and dashed 
curves. From left to right, each pair has fixed baryon charge 
B = 10 20 (red), 10 24 (black) and 10 33 (blue) respectively. The 
solid curves represent nuclear density cores while the dashed 
curves represent 100 times nuclear density. The light shaded 
(yellow) regions correspond to the nonrelativistic Boltzmann 
regime where the Boltzmann approximation (10) (dash-dotted 
(cyan) line) is valid. Only the B = 10 20 profiles visibly de- 
part from the one-dimensional approximation in this regime. 
The two upper curves in the top plot use the scale on the 
right and are the annihilation rates Yp s (right (green) curve) 
and r^ ir (left (blue) curve) normalized to the saturated value 
r?s(^ = oo ) = 4i^ 3 /(37rra 2 A: 2 ) (12). These curves comprise 
a range of cutoffs q ~ mot such the positronium annihilation 
rates (12) vary by 10%, and a range of incoming velocities 
10~ 3 c < v < 10~ 2 c. The scaling of the abscissa in the up- 
per figure is ln(z + z UR ) where z UR w 5 x 10 -11 so that the 
T = ultrarelativistic approximation (9) is linear (finely dot- 
ted (cyan) line). The abscissa are linearly spaced in the bottom 
figure. 
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so that the this approximation matches the numerical 
solution at = (raT/27r) 3/2 , (a slightly better approx- 
imation than used in [5]). 

All of our previous estimates about physical proper- 
ties of the antimatter nuggets - emission of radiation, 
temperature etc. - have been based on these approxima- 
tions using the region indicated in Fig. 1, it is clear that 
they work very well, as long as these regimes exist. The 
potential problem is that the nuggets could be highly 
charged or sufficiently small that the one-dimensional 
Boltzmann regime fails to exist with the density rapidly 
falling through the relevant density scales. We can see 
from the results in Fig. 1 that the nonrelativistic re- 
gion exists for all but the smallest nuggets: As long 
as B > 10 20 - as required by current detector con- 
straints [4] - then the one-dimensional approximation 
suffices to calculate the physical properties. 

The annihilation rates discussed in Sec. V have also 
been included in Fig. 1 to show that these effects only 
start once the densities have reached the atomic den- 
sity scale, which is higher than the Boltzmann regime. 
Thus these results are also insensitive to the size of the 
nugget and the domain-wall approximation suffices. It 
is clear that one must have a proper characterization of 
the entire density profile from nonrelativistic through 
to ultrarelativistic regimes in order to properly calculate 
the emissions: one cannot use simple analytic forms. 

We now discuss how charge equilibrium is estab- 
lished, and then apply our results to fix the relative nor- 
malization between the diffuse 1-20 MeV emissions 
and the 511 keV emissions from the core of our Galaxy. 
As we shall show, the relative normalization is now 
firmly rooted in conventional physics, and in agreement 
with the previous phenomenological estimate [6], pro- 
viding another validation of our theory for dark matter. 

C. Nugget Charge Equilibrium 

In order to determine the effective charge of the nuggets 
we must consider how equilibrium is obtained. For 
the matter nuggets, equilibrium is established through 
essentially static equilibrium with the surrounding 
ism plasma, however, for the antimatter nuggets, no 
such static equilibrium can be achieved. Instead, one 
must consider the dynamics of the following charge- 
exchange processes: 

1. Deposition of charge via interaction and annihilation of 
neutral ism components (primarily neutral hydrogen). 

2. Deposition of charge via interaction and annihilation of 
ionized ism components (electrons, protons, and ions). 

3. Evaporation of positrons from the antinugget's surface. 

First we consider impinging neutral atoms or 
molecules with velocity £> N/ i SM ~ 10 -3 c. Even if the 



antinugget is charged, these are neutral and will still 
penetrate into the electrosphere. Here the electrons 
will annihilate or be ionized leaving a positively charge 
nucleus with energy T ISM ~ 2^n^n,ism- As tne e ^ ec ~ 
trosphere consists of positrons, this charge cannot be 
screened, so the nucleus will accelerate to the core. At 
the core the nucleus either will penetrate and annihi- 
late, resulting in the diffused x-ray emissions discussed 
in [4] and providing the heat to fuel the microwave 
emissions [5], or will bounce off of the surface due to 
the sharp quark-matter interface. 9 




FIG. 2. Charge-exchange process (top diagram): A charged 
incoming proton (upper left) exchanges an up quark u for a 
down quark d with the antimatter nugget B, reflecting as a 
neutron (upper right). This is Zweig suppressed relative to 
the simple reflection process illustrated below, however, the 
neutron can escape from the system whereas the reflected pro- 
ton will be trapped by the electric fields. This will enhance the 
overall rate of charge-exchange reactions because the proton 
will continue to react again and again until either the charge 
exchange occurs, or it eventually annihilates. In order to ex- 
plain the relative intensities of the observed 511 keV and dif- 
fuse x-ray emissions from the core of our Galaxy [4], the ratio 
of the charge-exchange rate to the annihilation rate must be 
/ ~ 0.1-0.5. To calculate this from microscopic physics, how- 
ever, requires details of the nugget's surface and equation of 
state that are not yet known. Perhaps some model-insensitive 
estimates could be made, which would provide a highly non- 
trivial test of our theory based on microscopic physics. 



9 It is well known from quantum mechanics that any sufficiently 
sharp transition has a high probability of reflection of low energy 
particles. 
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Note that, if the nuggets are sufficiently charged such 
that kinetic energy of the incoming nucleus at the point 
of ionization is smaller than the electrostatic energy, 

T ISM ~ I^n^ism ~ eV < 6 -y ~ keV - 10 7 K, (11) 

then the charged nucleus will be unable to escape and 
will return to the core to either annihilate or undergo 
a charge-exchange reaction to become neutral, after 
which it may leave unimpeded by the electric field. For 
nuggets with B ~ 10 20 , 10 24 , and 10 33 respectively, this 
critical charge is Q/e ^> 1, 10 2 , and 10 4 respectively As 
shown in Table I, the charge established through evapo- 
ration greatly exceeds this critical charge: The charged 
nucleus will not be able to escape the antinugget unless 
it undergoes a charge-exchange reaction (see Fig. 2). 

For a single proton, such a charge-exchange reaction 
consists of an up quark being replaced by a down quark 
at the surface of the nugget, converting the proton to a 
neutron which can then escape. This process involves 
the exchange of a quark-antiquark pair and is a strong 
interaction process, but suppressed by the Zweig (ozi) 
rule. The overall ratio of charge-exchange to annihila- 
tion reactions is amplified by a finite probability of re- 
flection from the core boundary: this results in multiple 
bounces before annihilation. 

A better understanding of the details concerning the 
interaction between the proton and the quark-matter 
boundary is required to quantitatively predict the ratios 
of the rate of charge-exchange to the rate of annihilation 
and thus to confirm or rule out the suppression factor 
/ ~ 0.1-0.5 required to explain the relative strengths of 
the observed 511 keV and diffuse ~ 10 keV x-ray emis- 
sions seen from the core of the Galaxy [4]. This requires 
an estimate of the rates of reflection, annihilation, and 
charge-exchange; the elasticity of collisions; and the en- 
ergy loss through scattering. 

We emphasize that the corresponding calculations 
do not require any new physics: everything is rooted 
firmly in qed and qcd. They do, however, require solv- 
ing the many-body physics of the strong interaction in- 
cluding the incident nucleons and the quark antimatter 
interface (the phase of which may be quite complicated). 
A full calculation is thus very difficult, requiring insight 
into high-density qcd: estimates can probably be made 
using standard models of nuclear matter for the core, 
but are beyond the scope of the present paper. 

In any case, the nucleus will certainly deposit its 
charge on the antinugget, and one may neglect the in- 
teractions with neutral ism components for the purpose 
of establishing charge neutrality. 

Instead, we must consider the charged components. 
Using the same argument (11) one can see that the 
Coulomb barrier of the charged nuggets will be high 
enough to prevent electrons from reaching the electro- 
sphere in all but the very hot ionized medium (vhim), 



which occupies only a small fraction of the ism in the 
core (see [47] for example). Protons however, will be 
able to reach the core to annihilate with a cross sec- 
tion rsj 4nR 2 (this is not substantially affected by the 
charge). Thus, positive charge can be deposited at a 
rate ~ 47tR 2 v N + /lSM n N + 

,ism where n N + ISM is the density 
of the ionized components in the ism, which is typically 
10 _1 to 10 -2 of the total density n ISM ~ 1 cm -3 

This rate of positive charge accumulation must match 
the evaporation rate of the positrons from the elec- 
trosphere (see Eq. (7)), giving the boundary condi- 
tions (8b), and the resulting charges summarized in Ta- 
ble I. (These are consistent with the estimates in [18] 
and the universal upper bound [48].) 

Note that the rate of evaporation is much less than 
the rate of proton annihilation and carries only ~ 1 eV 
of every per particle compared to the nearly 2 GeV of 
energy deposited. Thus, evaporation does not signifi- 
cantly cool the nuggets: the thermal radiation from the 
Boltzmann regime discussed in [5] dominates. 



V. DIFFUSE GALACTIC EMISSIONS 

Previous discussions of quark (anti)nugget dark matter 
considered the electrosphere in only two limiting cases: 
the inner high-density ultrarelativistic regime and the 
outer low-density Boltzmann regime. While this anal- 
ysis was sufficient to qualitatively discuss the different 
components of the spectrum, it did not allow for their 
comparison in any level of detail. In particular, com- 
paring the relative strengths of the 511 keV line emis- 
sion (emitted entirely from the nonrelativistic regime) 
with that of the MeV continuum, (emitted from the 
relativistic regime) required introducing a phenomeno- 
logical parameter % « 0.1 [6] to express the relative 
rates of direct annihilation to positronium formation. 
This parameter is sensitive to the density profile at all 
scales. Using the detailed numerical solutions of the 
previous section we can directly compute the relative 
intensities of the two emission mechanisms. Our calcu- 
lation demonstrates that this value of % is supported by 
a purely microscopic calculation rooted firmly in qed 
and well-understood many-body physics. The agree- 
ment between the calculated value of % and the phe- 
nomenological value required to fit the observations 
provides another important and nontrivial test of our 
dark-matter proposal. In addition, using the same nu- 
merical solutions of the previous section, we compute 
the spectrum in the few MeV region which could not be 
calculated in [6] with only the ultrarelativistic approxi- 
mation for the density profile. 
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A. Observations 

1. The 511 keV Line 

Spi/ integral has detected a strong 511 keV signal from 
the Galactic bulge [49]. A spectral analysis shows this 
to be consistent with low-momentum e + e~ annihila- 
tion through a positronium intermediate state: About 
one quarter of positronium decays emit two 511 keV 
photons, while the remaining three-quarters will de- 
cay to a three photon state producing a continuum 
below 511 keV. Both emissions have been seen by 
spi/ integral with the predicted ratios. 

The 511 keV line is strongly correlated with the Galac- 
tic centre with roughly 80% of the observed flux coming 
from a circle of half angle 6°. There also appears to be a 
small asymmetry in the distribution oriented along the 
Galactic disk [50]. 

The measured flux from the Galactic bulge is found 
to be dO/dQ ~ 0.025 photons cm -2 s _1 sr" 1 [51]. After 
accounting for all known Galactic positronium sources 
the 511 keV line seems too strong to be explained by 
standard astrophysical processes. These processes, as 
we understand them, seem incapable of producing a 
sufficient number of low-momentum positrons. Several 
previous attempts have been made to account for this 
positron excess. Suggestions have included both modi- 
fications to the understood spectra of astrophysical ob- 
jects or positrons which arise as a final state of some 
form of dark-matter annihilation. At this time there is 
no conclusive evidence for any of these proposals. 

If we associate the observed 511 keV line with the 
annihilation of low momentum positrons in the elec- 
trosphere of an antiquark nugget, these properties are 
naturally explained. The strong peak at the Galactic 
centre and extension into the disk must arise because 
the intensity follows the distribution py(r)p BM (r) of vis- 
ible and dark-matter densities. This profile is unique to 
dark-matter models in which the observed emission is 
due to matter-dark matter interactions and should be 
contrasted with the smoother ~ Pdm(^) profile for pro- 
posals based on self-annihilating dark-matter particles, 
or the rsj p DM (r) profile for decaying dark-matter pro- 
posals. The distribution pv(r)p BM (r) obviously implies 
that the predicted emission will be asymmetric, with 
extension into the disk from the Galactic center as it 
tracks the visible matter. There appears to be evidence 
for an asymmetry of this form [50]. In our proposal, 
no synchrotron emission will occur as the positrons are 
simply an integral part of the nuggets rather than be- 
ing produced at high energies with only the low-energy 
components exposed for emission. Contrast this with 
many other dark matter based proposals where the rel- 
atively high-energy positrons produced from decaying 
or annihilating dark-matter particles produce strong 



synchrotron emission. These emissions are typically in 
conflict with the strong observational constraints [21]. 



2. Diffuse MeV scale emission 

The other component of the Galactic spectrum we dis- 
cuss here is the diffuse continuum emission in the 1 - 
30 MeV range observed by comptelcrgo. The interpre- 
tation of the spectrum in this range is more complicated 
than that of the 511 keV line as several different astro- 
physical processes contribute. 

The conventional explanation of these diffuse emis- 
sions is that of gamma rays produced by the scattering 
of cosmic rays off of the interstellar medium, and while 
detailed studies of cosmic ray processes provide a good 
fit to the observational data over a wide energy range 
(from 20 MeV up to 100 GeV), the predicted spectrum 
falls short of observations by roughly a factor of 2 in the 
1-20 MeV range [52]. 

Background subtraction is difficult, however, and es- 
pecially obscures the spatial distribution of the MeV ex- 
cess. However, the excess seems to be confined to the 
inner Galaxy (/ = 330° - 30°, \b\ = 0° - 5°), [52] with a 
negligible excess from outside of the Galactic centre. 



3. Comparison 

As our model predicts both of these components to have 
a common source, a comparison of these emissions pro- 
vides a stringent tests of the theory. In particular, the 
morphologies, spectra, and relative intensities of both 
emissions must be strongly related. 

The inferred spatial distribution of these emissions 
is consistent: both are concentrated in the core of the 
Galaxy. Unfortunately, this is not a very stringent test 
due to the poor spatial resolutions of the present ob- 
servations. The prediction remains firm: If the mor- 
phology of the observations can be improved, then a 
full subtraction of known astrophysical sources should 
yield a diffuse MeV continuum with spatial morphol- 
ogy identical to that of the 511 keV. 

Present observations, however, do allow us to test the 
intensity and spectrum predicted by the quark nugget 
dark-matter model. The model predicts the intensity 
to be proportional to the 511 keV flux with calculable 
coefficient of proportionality. Therefore, the integral 
data may be used to fix the total diffuse emission flux, 
removing the uncertainties associated with the line of 
sight averaging which is the same for both emissions. 

The resulting spectrum and intensity were previously 
discussed in [24] and [6], but the 511 keV line emission 
was estimated using the low-density Boltzmann approx- 
imation, while the MeV emissions were estimated using 
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the high-density ultrarelativistic approximation. The 
proportionality factor linking these two was treated as a 
phenomenological parameter, %, which required a value 
of ^ ~ 0.1 [6] in order to explain the observations. Here, 
equipped with a complete density profile, we calculate 
the value of this parameter from first principles show- 
ing that it is indeed consistent with the observations, 
providing yet another highly nontrivial verification of 
the quark antimatter nugget dark-matter proposal. 



B. Annihilation rates 

As incident electrons enter the electrosphere of 
positrons, the dominant annihilation process is through 
positronium formation with the subsequent annihila- 
tion producing the 511 keV emission. As the remaining 
electrons penetrate more deeply, they will encounter a 
higher density of more energetic positrons: Direct anni- 
hilation eventually becomes the dominant process. The 
maximum photon energy, ~ 10 MeV, is determined by 
the Fermi energy of the electrosphere at the deepest 
depth of electron penetration. We stress that this scale is 
not introduced in order to explain the comptel/ crgo 
gamma-ray excess: Our model necessarily produces a 
strong emission signature at precisely this energy. 

In the following section we integrate the rates for 
these two processes over the complete density profile, 
allowing us to predict the relative intensities and spec- 
tral properties of the emissions. 



1. Positronium Formation 

In principle, we only need to know the e + e~ annihila- 
tion cross section at all centre of mass momenta. At 
high energies, this is perturbative and one can use the 
standard qed result for direct e + e~ —> 2j emission. At 
low energies, however, one encounters a strong reso- 
nance due to the presence of a positronium state the 
greatly enhances the emission. This resonance renders 
a perturbative treatment invalid and must be dealt with 
specially, thus we consider these as two separate pro- 
cesses. Our presentation here will be abbreviated: de- 
tails may be found in [24]. 

While the exact positronium formation rate - 
summed over all excited states - is not well established, 
it is clear that the rate will fall off rapidly as the centre of 
mass momentum moves away from resonance. A sim- 
ple estimate suggests that for momenta p > moc the for- 
mation rate falls as ~ p -4 . To estimate the rate, we thus 
make a cutoff at q « moc as the upper limit for positro- 
nium formation: for large center of mass momenta we 
use the perturbative direct-annihilation approximation. 
The scale here is set by the Bohr radius = 2(ma)~ l 



for the positronium bound state: If a low-momentum 
e + e~ pair pass within this distance, then the probability 
of forming a bound state becomes large, with a natural 
cross section ap s ~ it a 2 . The corresponding rate is 




where n(p) is the momentum distribution of the 
positrons in the electrosphere, and v ~ oc is the incom- 
ing electron velocity. The second expression represents 
the two limits (valid at low temperature): a) of low den- 
sity when the Fermi momentum pp < q and all states 
participate, resulting in a factor of the total density n, 
and b) of high density where the integral is saturated 
by q 3 /3n 2 « (moc) 3 /3rc 2 . 

As discussed in Sec. IV C, and in more detail in Ap- 
pendix B, electrons will be able to penetrate the charged 
antimatter nuggets in spite of the strong electric field. 
Initially the electrons are bound in neutral atoms. Once 
ionized, the density is sufficiently high that the charge 
is efficiently screened with a Debye screening length Ap 
that is much smaller than the typical de Broglie wave- 
length A = h/p of electrons. Thus, the electric fields 

- although quite strong in the nugget's electrosphere 

- will not appreciably effect the motion. The binding 
of electrons in neutral atoms complicates the analysis 
slightly, but the binding energy - on the eV scale - will 
not significantly alter the qualitative nature of our esti- 
mates. For a precision test of the emission properties, 
this will need to be accounted for. Here we include 
bands comprising ±10% relative variation in the over- 
all positronium annihilation rate to show the sensitivity 
to this uncertainty. 

2. Direct-annihilation Rates 

While positronium formation is strongly favoured at 
low densities due to its resonance nature, deep within 
the electrosphere the rapidly growing density of states 
at large momentum values result in a cross section 
characterized by the perturbative direct-annihilation 
process. Conceptually one can imagine an incident 
Galactic electron first moving through a Fermi gas of 
positrons with roughly atomic density with a relatively 
large probability of annihilation through positronium 
to 511 keV photons. A small fraction survive to pen- 
etrate to the inner highOdensity region where direct 
annihilation dominates. The surviving electrons then 
annihilate with a high-energy positron near to the in- 
ner quark-matter surface releasing high-energy pho- 
tons. The spectrum of these annihilation events will 
be a broad continuum with an upper cutoff at an en- 
ergy scale set by the Fermi energy of the positron gas at 
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Distance from nugget surface: z = r — R (in cm) 

FIG. 3. Electron survival fraction, n e (r)/noo (16), of an in- 
coming electron with velocities v = 0.01c [leftmost light gray 
(red) band], v = 0.005c [middle (green) band] and v = 0.001c 
[rightmost dark gray (blue) band] from left to right, respec- 
tively. The thickness of the bands includes a ±10% variation 
in the positronium annihilation rate (12). The local Fermi mo- 
mentum pp is shown along the top and with vertical dotted 
lines including the cutoff scale q w 3.7 keV from (12). The yel- 
low shaded Boltzmann region and the abscissa scaling are the 
same as in the top of Fig. 1. Note that annihilation happens 
well within the electrosphere, so the finite size of the nugget 
is irrelevant. 



the maximum penetration depth of the electrons. The 
spectral density for direct e + e~ annihilation at a given 
chemical potential was calculated in [6]: 
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where Ep = ^ p 2 + ra 2 is the energy of the positron in 
the rest frame of the incident (slow-moving) electron 
and cv is the energy of the produced photons. The anni- 
hilation rate at a given density, r^ ir (n^), is obtained by 
integrating over allowed final state photon momentum. 
This was previously done in the T —> limit where the 
integral may be evaluated analytically, at nonzero tem- 
peratures the rate must be evaluated numerically. 



C. Spectrum and Branching Fraction 

To determine the full annihilation spectrum we first de- 
termine the fraction of incident electrons that can pene- 



trate to a given radius r in the electrosphere (see Fig. 3. 
We then integrate the emissions over all regions. 

Consider an incident beam of electrons with density 
ftoo and velocity v. As they enter the electrosphere of 
positrons, the electrons will annihilate. The survival 
fraction n e (r)/ftoo will thus decrease with a rate pro- 
portional to r(r) = Tp s (r) + T dir (r) which depends cru- 
cially on the local density profile n[ji(r)] calculated in 
Sec. IV B: 

= v = -Ttfv^neir). (15) 

Integrating (15), we obtain the survival fraction: 



n e (r 



exp 



T(r) 



(16) 



This is shown in Fig. 3. One can clearly see that in the 
outer electrosphere, positronium formation - indepen- 
dent of v - dominates the annihilation. Once the density 
is sufficiently high (pp > 1 MeV), the direct-annihilation 
process dominates, introducing a dependence on the ve- 
locity v of the incident particle. 

The initial velocity v ~ 10 _3 c is determined by the 
local relative velocity of the nuggets with the surround- 
ing ism. This will depend on the temperature of the 
ism, but the positronium annihilation rate is insensitive 
to this. As we mentioned previously, most of the elec- 
trons in the interstellar medium are bound in neutral 
atoms, either as neutral hydrogen HI, or in molecular 
form H2. (The ionized hydrogen HII represents a very 
small mass fraction of interstellar medium.) These neu- 
tral atoms and molecules will have no difficulty enter- 
ing the electrosphere. 

The remaining bound electrons that do not annihilate 
through positronium will ionize once they reach denser 
regions, and will acquire a new velocity set by a combi- 
nation of the initial velocity v ~ 10 _3 c and the atomic 
velocity v ~ oc ~ 10~ 2 c imparted to the electrons as 
they are ionized from the neutral atoms: The latter will 
typically dominate the velocity scale. This only occurs 
in sufficiently dense regions where the Debye screen- 
ing discussed in Appendix B becomes efficient. Hence, 
the electric fields will not significantly alter the motion 
of the electrons after ionization. The direct-annihilation 
process depends on the final velocity v; as the dominant 
contribution comes from ionization, this will remain rel- 
atively insensitive to the ism. During ionization, some 
fraction (roughly half) of the electrons will move away 
from the core, but a significant portion will travel with 
this velocity v toward the denser regions. 

Two other features of Fig. 3 should be noted. First 
is the value of the survival fraction % ~ 0.1 at which 
direct annihilation dominates [see Eq. 19]. This is the 
value that was postulated phenomenologically in [6] in 
order to explain the relative intensities of the 511 keV 
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and direct-annihilation spectra. Here this value results 
naturally from a microscopic calculation in a highly 
nontrivial manner that depends on the structure of 
the density profile at all scales up to p f = 1 MeV. 
Second is the rapid increase in density near the core 
quickly extinguishes any remaining electrons. As such, 
even if the chemical potential at the surface is larger - 
jiR ~ 100 MeV or so - virtually no Galactic electrons 
penetrate deeply enough to annihilate at these energies. 
For nonrelativistic electrons the maximum energy scale 
for emission is quite generally set at the ~ 20 MeV scale, 
depending slightly on v. 

Having established the survival fraction as a function 
of height it is now possible to work out the spectral 
density that will arise from the e + e~ annihilations. At 
a given height we can express the number of annihila- 
tions through a particular channel as, 

^=n e (r)v- 1 T(r). (17) 

Integrating this expression over all heights will then 
give the total fraction of annihilation events proceeding 
via positronium f 5 n, and via direct annihilation fyieV- 
(The numerical values are given for q = moc and are 
quite insensitive to v.) 

f 5n = r dr ?M v -i Ips ( r ) « o.9, (18a) 

JR ft co 

fuev = r dr ^v-'TaUr) « 0.1. (18b) 

JR ft co 

Note that in both cases, the overall normalization de- 
pends only on the total rate of electron collisions with 
dark matter through vrioo. As such, the ratio of 511 keV 
photons to MeV continuum emission is independent of 
the relative densities (though it will show some depen- 
dence on the local electron velocity distribution). Nu- 
merically, with a ±10% variation in the positronium an- 
nihilation rate, we find 

x = f^L « 0.05 - 0.2, (19) 

7511 

which is quite insensitive to the velocity v. This ratio 
was introduced as a purely phenomenological param- 
eter in [6] to explain the observations. Here we have 
calculated from purely microscopic considerations that, 
for a wide range of nugget parameters, the required 
value x ~ 0-1 [6] arises quite naturally. 

We now have everything needed to compute the spec- 
tral density of the MeV continuum. The spectral density 
at a fixed chemical potential is given by Eq. (13): this 
must now be averaged along the trajectory of the in- 
coming electron weighted by the survival fraction (16) 
and the time spent in the given region of the trajectory 
(set by the inverse velocity v~ l )\ 

ductal _ r df? -i Mr)dlW) (2Q) 
dco Jr ft co dco dt 




uj (MeV) 



FIG. 4. Spectral density (scaled by co 2 to compare with [52] 
of photons emitted by an electron annihilating on antiquark 
nuggets with incoming velocities v = 0.01c [uppermost (red) 
band], v = 0.005c [middle (green) band] and v = 0.001c [low- 
est (blue) band to the lower left] from right to left, respec- 
tively, including the cosmic ray background determined in 
[52] (dotted line). The thickness of the bands includes a ±10% 
variation in the positronium annihilation rate (12). The over- 
all normalization is fixed to the observationally unrelated 511 
keV line as discussed below (21). The three error bars are 
the comptel data points. Note: This spectrum should still be 
interpreted as a qualitative effect - a detailed calculation of 
the ionization and hence distribution of the velocity v must 
be performed to yield a quantitative prediction. The general 
structure and magnitude, however, can be trusted as these de- 
pend on the overall density profile which we have carefully 
modeled. (Compare with Fig. 5, for example, which uses only 
the ultrarelativistic density profile: The resulting intensity is 
2 orders of magnitude too large.) 

The resulting spectrum is shown in Fig. 4 and is sensi- 
tive to both the incoming velocity v (which will depend 
on the local environment of the antimatter nugget) and 
the overall normalization of the positronium annihila- 
tion rate (12). (The latter is fixed in principle, but re- 
quires a difficult in-medium calculation to determine 
precisely.) These two parameters are rather orthogonal. 
The velocity v determines the maximum depth of pene- 
tration, and hence the maximum energy of the emitted 
photons (as set by the highest chemical potential at the 
annihilation point): If the electron velocity is relatively 
low (v < 100 km/ s « 0.0003c), almost all annihilations 
happen immediately and the MeV continuum will fall 
rapidly beyond 5 MeV. As the velocity increases the elec- 
trons are able to penetrate deeper toward the quark sur- 
face and annihilate with larger energies. In contrast, the 
details of the positronium annihilation do not alter the 
spectral shape, but do alter the overall normalization. 

As already mentioned earlier, the parameter % ~ 
0.1 is calculated here from purely microscopic physics. 
Therefore, it provides highly nontrivial verification of 
the entire proposal. Also, the profile function from the 
previous section is computed at all scales, allowing us 
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to calculate the photon spectrum down to, and below, 
the electron mass. This calculation could not be per- 
formed in the previous analysis [6] with only the ultra- 
relativistic expression. 

We stress here that the result shown in Fig. 4 ex- 
hibits several nontrivial features arising from very gen- 
eral characteristics of the proposed emission model. 
The wide range of positron energies within the electro- 
sphere implies a broad emission profile with a width 
Aco ~ 5-10 MeV. The initial growth of emission 
strength with to is due to the increasing density of states 
as a function of depth. Above ~ 10 MeV, the emission 
becomes suppressed by the inability of Galactic elec- 
trons to penetrate to depths where the positrons have 
this energy. While the exact details may vary with pa- 
rameters, such as the local electron velocity and precise 
rate of positronium formation, the general spectral fea- 
tures are inescapable consequences of our model, allow- 
ing it to be tested by future, more precise, observations. 



D. Normalization to the 511 keV line 

To obtain the observed spectrum, one needs to aver- 
age (20) along the line of site over the varying matter 
and dark-matter density and velocity distributions. Nei- 
ther of these is known very well, so to check whether or 
not the prediction is significant, we fix the average rate 
of electron annihilation with the observed 511 keV line 
which our model predicts to be produced by the same 
process. As the intensity of the 511 keV line (resulting 
from the two-photon decay of positronium in the 1 So 
state) has been measured by spi/integral along the 
line of sight toward the core of the Galaxy (see [53] for 
a review), we can use this to fix the normalization of the 
MeV spectrum along the same line of sight. The total 
intensity is given in terms of Eq. (18a): 



dO _ C, 
dQ - 4 /5n 



0.025 



photons 
cm 2 s sr 



(21) 



where the factor of 4 accounts for the three-quarters 
of the annihilation events that decay via the 3 Si chan- 
nel (also measured, but not included in the line emis- 
sion). This fixes the normalization constant C « 
0.11 events cm -2 s _1 sr -1 . The predicted contribution 
to the MeV continuum must have the same morphol- 
ogy and, thus, an integrated intensity along the same 
line of sight must have the same normalization factor. 
This normalization has been used in Fig. 4 to compare 
the predicted spectrum with the unaccounted for excess 
emission detected by comptel ([52]). 

The exact shape of the spectrum will be an average 
of the components shown over the velocity distribution 
of the incident matter. Our process of normalization is 
unable to remove this ambiguity because the predicted 



511 keV spectral properties are insensitive to this. It is 
evident, however, that MeV emission from dark-matter 
nuggets could easily provide a substantial contribution 
to the observed MeV excess. Note also, that the ex- 
cess emission can extend only to 20 MeV or so. This is 
completely consistent with the more sophisticated back- 
ground estimates discussed in [54] which can fit almost 
all aspects of the observed spectrum except the excess 
between 1 and 30 MeV predicted by our proposal. 



VI. CONCLUSION 

Solving the relativistic Thomas-Fermi equations, we de- 
termined the charge and structure of the positron elec- 
trosphere of quark antimatter nuggets that we postu- 
late could comprise the missing dark-matter in our Uni- 
verse. We found the structure of the electrosphere to 
be insensitive to the size of the nuggets, as long as they 
are large enough to be consistent with current terrestrial 
based detector limits, and hence, can make unambigu- 
ous predictions about electron annihilation processes. 

To test the dark-matter postulate further, we used 
the structure of this electrosphere to calculate the an- 
nihilation spectrum for incident electronic matter. The 
model predicts two distinct components: a 511 keV 
emission line from decay through a positronium inter- 
mediate and an MeV continuum emission from direct- 
annihilation processes deep within the electrosphere. 
By fixing the general normalization to the measured 511 
keV line intensity seen from the core of the Galaxy, our 
model makes a definite prediction about the intensity 
and spectrum of the MeV continuum spectrum without 
any additional adjustable model parameters: Our pre- 
dictions are based on well-established physics. 

As discussed in [30], a difficulty with most other 
dark-matter explanations for the 511 keV emission is 
to explain the large ^100% observed positronium frac- 
tion - positrons produced in hot regions of the Galaxy 
would produce a much smaller fraction. Our model 
naturally predicts this observed ratio everywhere. 

A priori, there is no reason to expect that the predicted 
MeV spectrum should correspond to observations: typi- 
cally two uncorrelated emissions are separated by many 
orders of magnitude. We find that the phenomenolog- 
ical parameter % (19) required to explain the relative 
normalization of MeV emissions arises naturally from 
our microscopic calculation. This is highly nontriv- 
ial because it requires a delicate balance between the 
two annihilation processes from the semirelativistic re- 
gion of densities that is sensitive to the semirelativis- 
tic self -consistent structure of the electrosphere outside 
the range of validity of the analytic ultrarelativistic and 
nonrelativistic regimes. (See Fig. 1 and 3). 

If the predicted emission were several orders of mag- 



15 



nitude too large, the observations would have ruled 
out our proposal. If the predicted emissions were too 
small, the proposal would not have been ruled out, but 
would have been much less interesting. Instead, we 
are left with the intriguing possibility that both the 511 
keV spectrum and much of the MeV continuum emis- 
sion arise from the annihilation of electrons on dark- 
antimatter nuggets. While not a smoking gun - at least 
until the density and velocity distributions of matter 
and dark-matter are much better understood - this pro- 
vides another highly nontrivial test of the proposal that, 
a priori, could have ruled it out. 

Both the formal calculations and the resulting struc- 
ture presented here - spanning density regimes from ul- 
trarelativistic to nonrelativistic - are similar to those rel- 
evant to electrospheres surrounding strange-quark stars 
should they exist. Therefore, our results may prove 
useful for studying quark star physics. In particular, 
problems such as bremsstrahlung emission from quark 
stars originally analyzed in [55] (and corrected in [56]) 
that uses only ultrarelativistic profile functions. The re- 
sults of this work can be used to generalize the corre- 
sponding analysis for the entire range of allowed tem- 
peratures and chemical potentials. Another problem 
which can be analyzed using the results of the present 
work is the study of the emission of energetic electrons 
produced from the interior of quark stars. As advo- 
cated in [57], these electrons may be responsible for 
neutron star kicks, helical and toroidal magnetic fields, 
and other important properties that are observed in a 
number of pulsars, but are presently unexplained. 

Finally, we would like to emphasize that this mech- 
anism demonstrates that dark matter may arise from 
within the standard model at the qcd scale, 10 and that 
exotic new physics is not required. Indeed, this is nat- 
urally suggested by the "cosmic coincidence" of almost 
equal amounts of dark and visible contributions to the 
total density Otot = 1.011(12) of our Universe [16]: 



^dark-energy • Q 



dark-matter 



visible 



17:5:1. 



The dominant baryon contribution to the visible portion 
^visible ~ has an obvious relation to qcd through 
the nucleon mass ra N oc A QCD (the actual quark masses 
arising from the Higgs mechanism contribute only a 
small fraction to ra N ). Thus, a qcd origin for the dark 
components would provide a natural solution to the 
extraordinary "fine-tuning" problem typically required 
by exotic high-energy physics proposals. Our proposal 
here solves the matter portion of this coincidence. For 
a proposal addressing the energy coincidence we refer 
the reader to [58] and references therein. 11 



10 The axion is another dark-matter candidate arising from the qcd 
scale, but with fewer observational consequences. 

11 The idea concerns the anomaly that solves the famous axial U(1)a 
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end of appendix 1) where it was emphasized that [45] 
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Appendix A: Density Functional Theory 

We start with a Density Functional Theory (dft) formu- 
lated in terms of the thermodynamic potential: 

O = £2/i / d 3 f ipt(?)(e_ mt - p)ip i (r) + 

i J 

J d 3 rV ext (r)n(r) + J d 3 re xc n(r) n(r) + 

J^r^^p^-TEfM, (Ala) 



+ 



Arte 2 



where e v = \/p 2 + ra 2 is the relativistic energy of the 
electrons, £ X c(^) is the exchange energy. The density is 



n(r)=2Y J fif i (r)Ur) l 



(Alb) 



where we have explicitly included the spin degeneracy, 
and used relativistic units where h = c = 1 and e 2 = oc. 

We vary the potential with respect to the occupation 
numbers /; and the wave functions ipj subject to the con- 
straints that the wave functions be normalized, 

J d 3 rxp}(r)xp j (r)=S ij , (Ale) 



problem, giving rise to an rj r mass that remains finite, even in the 
chiral limit. Under some plausible and testable assumptions about 
the topology of our Universe, the anomaly demands that the cosmo- 
logical vacuum energy depend on the Hubble constant H and qcd 
parameters as p VE ~ Hm q {qq} fm^i « (4 x 10~ 3 eV) 4 - tantalisingly 
close to the value p DE = [1.8(1) x 10~ 3 eV] 4 observed today [16]. 
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and that the occupation numbers satisfy Fermi-Dirac 
statistics. 12 This yields the Kohn-Sham equations for 
the electronic wave functions: 

e (-ihV) ~V + ^eff (»*)]#(**) = 

f.= 1 

* 1+e £,/T- 

where (we take e to be positive here) 

V e& (f) = e0(f ) + Vext(r) + e«[n(r)] + n(f)e' xc [n(f)] 
and 

is the electrostatic potential obeying Poisson's equation, 
V 2 (p(r) —^Tcen(r). 

1. Thomas-Fermi Approximation. 

If the effective potential V e ft(r) varies sufficiently 
slowly, then it is a good approximation to replace it 
locally with a constant potential. The Kohn-Sham equa- 
tions thus become diagonal in momentum space, i = fc, 
ipk oc e lkr and may be explicitly solved: 



n[r 



d 3 k 



~H + Vef f (r), 
1 1 



(2zr) 3 Vl + eW /T l + e - E ^ f ) /T 



where we have also explicitly included the contribu- 
tion from the antiparticles. We piece these homoge- 
neous solutions together at each point r and find the 
self-consistent solution that satisfies Poisson's equation. 
This approximation is a relativistic generalization of the 
Thomas-Fermi approximation. 

In principle, the dft method is exact [59], however, 
the correct form for £ xc is not known and could be 
extremely complicated. Various successful approxima- 
tions exist but for our purposes we may simply ne- 
glect this. The weak electromagnetic coupling constant 
oc ~ 1/137 and Pauli exclusion principle keep the elec- 
trons sufficiently dilute that the many-body correlation 
effects can be neglected until the density is high in the 
sense that oc 2 n ~ oc 2 fi 3 /3n 3 > m 3 , which corresponds to 
]i > 50 MeV. 

Near the nugget core, many-body correlations may 
become quantitatively important. For example, the ef- 
fective mass of the electrons in the gas is increased 



12 This is most easily realized by introducing the single-body density 
matrix p — 1 — C p T C where C is the charge-conjugation matrix. 



by about 20% when ]i « 25 MeV and doubles when 
]i « 100 MeV (see, for example, [60]). These effects, 
however, will not change the qualitative structure, and 
can be quite easily taken into account if higher accuracy 
is required. We also note that, formally, the Thomas- 
Fermi approximation is only valid for sufficiently high 
densities. However, it gives correct energies within fac- 
tors of order unity for small nuclei [44] and is known to 
work substantially better for large nuclei [61]. As long 
as we do not attempt to use it in the extremely low- 
density tails, it should give an accurate description. 



2. Analytic Solutions 

There are several analytic solutions available if we con- 
sider the one-dimensional approximation, neglecting 
the curvature term 2ji' (r) /r in (5), which is valid close 
to the nugget where z = r — R <C R, the distance from 
the nugget core, is less than the radius of the nugget. 



a. Ultrarelativistic Regime 

The first is the ultrarelativistic approximation where ]i 
and /or T are much larger than m and the limit m —¥ 
can be taken. In this case, we may explicitly evaluate 
(see also [40-42]), 13 



3tt 2 ^ 3 " 



(A5) 



In the domain-wall approximation, this admits an exact 
solution with a typical length scale of z UR : 



]i{z f T) 



ZuR 




(A6) 



(A7) 



13 In this limit, then integrals have a closed form: 

n 2 Jo 

3 

dx 



l + e (p-p)/T l+ e {p+}i)/T 



tt2 



f 

Jo 



l + e x-}i/T l + e x+}i/l 



_ T 3 



where Li s (z) is the Poly logarithm 

oo 



r(3) [Li 3 (-e' l/T ) - Li 3 (-e-'' /T )l (A3) 



(A4) 
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In our case, X ~ 1 eV « m so we can also take the 
T -> limit to obtain (9): 



Wur(z) 



4 



3tt2(1+z/z ur )3' 




This solution persists until ]i ^ m, which occurs at a 
distance 



v m 



fr. Boltzmann Regime 

Once the chemical potential is small enough that 
e ji/T ^ e m//T , we may neglect the degeneracy in the 
system, and write 



n Mlexp(^/T)<exp(m/T) ~ n^ IT . 

This occurs for a density of about 



2n 



3/2 



(A10) 



(All) 



and lower. If the one-dimensional approximation is still 
valid, then another analytic solution may be found: 



n z 



2a: — 



(A12) 



from which we obtained (10). The shift Zb must be de- 
termined from the numerical profile at the point where 
n ( z B) = n B- Note that z a <C z B , so this regime is valid 
and persists until z ~ 0.1K, at which point the one- 
dimensional approximation breaks down. It is in this 
Boltzmann regime where most of the important radia- 
tive processes take place [5]. 

3. Numerical Solutions 

The main technical challenge in finding the numerical 
solution is to deal effectively with the large range of 
scales: T ~ 1 eV <C m ~ 500 keV < ji R ~ 25 MeV. For 
example, the density distribution (4b) is prone to round- 
off error, but the integral can easily be rearranged to 
give a form that is manifestly positive: 



1 f°° 



p 2 sinh (^) 



( ^±^)+cosh(f 



-. (A13) 



The derivative may also be safely computed. Let 



>/p 2 + m 1 



T 



(A14) 



to simplify the expressions. These can each be com- 
puted without any round-off error. The first derivative 
presents no further difficulties: 



1 y~ 2 1 



cosh A cosh B 



(cosh A + coshB) 



(A15) 



The differential equation is numerically simplified if we 
change variables to logarithmic quantities. We would 
also like to capture the relevant physical characteristics 
of the solution, known from the asymptotic regimes. 
Close to the nugget, we have ]i — }Ir/(1 + z/zq), we 
introduce an abscissa logarithmic in the denominator 



(A16) 



The dependent variable should be logarithmic in the 
chemical potential, so we introduce b = — ln(^/^o)- We 
thus introduce the following change of variables: 



r = £ + z (e fl -l), 
r-R 



(A17a) 



In 1 



ln^-. (A17b) 



With the appropriate choice of scales Zq describing the 
typical length scale at the wall r = R and fiQ ~ these 
form quite a smooth parametrization. The resulting sys- 
tem is 



&(«) 



47r az 2e 2fl ^l. (A18) 



It is imperative to include a full numerical solution 
to the profile in order to obtain the proper emission 
spectrum. Using only the ultrarelativistic approxima- 
tion (A8) produces a spectrum (Fig. 5) 2 orders of mag- 
nitude too large, in direct contradiction with the obser- 
vations [52]. The actual prediction depends sensitively 
on a subtle - but completely model-independent - bal- 
ance between the ultrarelativistic, relativistic, and non- 
relativistic regimes. The consistency between the pre- 
dicted spectrum shown in Fig 4 and the observations is 
a highly nontrivial test of the theory. 



Appendix B: Debye Screening In the Electrosphere 

Here we briefly discuss the plasma properties inside 
of the nugget's electrosphere. The main point is that 
the Debye screening length Ap is much smaller than 
the typical de Broglie wavelength A = ti/p of elec- 
trons. Thus, the electric fields - although quite strong in 
the nugget's electrosphere - will not appreciably effect 
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FIG. 5. Incorrect spectral density obtained by using a purely 
ultrarelativistic approximation for the profile (A8). Note that 
this is two orders of magnitude larger than the spectrum ob- 
tained from the full numerical solution (Fig. 4) and, if cor- 
rect, would have easily ruled out our proposal. This serves 
to demonstrate the highly nontrivial nature of the predicted 
spectrum shown in Fig. 4 that is consistent with the observa- 
tions [52]. 



the motion of electrons within the electrosphere as the 
charge is almost completely screened on a scale ~ Ad- 

This screening effect was completely neglected in [45], 
which led the authors to erroneously conclude that 
all electrons will be repelled before direct annihilation 
can proceed. The overall charge (see Table I) will re- 
pel incident electrons at long distances as discussed in 
Sec. IV C, but neutral hydrogen will easily penetrate the 
electrosphere, at which point the screening becomes ef- 
fective, allowing the electrons to penetrate deeply and 
annihilate as discussed in Sec. V. 

To estimate the screening of an electron, we solve the 
Poisson equation with a S 3 (r — Fq) function describing 
the electron at position tq (see also 2), 



-4rce 



n(r) — 5 (r — vq) 



(Bl) 



where (p(r) is the electrostatic potential and n(r) is 
the density of positrons. As before, we exchange (p 
for ]i, the effective chemical potential (3). Now let 
no(r) and }io(r) be the solutions discussed in Sec. IV B 
without the potential S 3 (r — vq). We may then de- 



scribe the screening cloud by ji(r) 
n(r) = fto(r) + 5 n {r) where 



po(r) +Su{r) and 



V 2 6 u (r) 



4zre z 



S n (r) -S 3 (r-r 
dn[p] 



Area 



dp. 



(B2) 



where n[p] is given by (4b). If the density is sufficiently 
large compared to the screening cloud deviations, we 



may take dn/dp to be a constant, in which case we may 
solve (B2) analytically with the boundary conditions: 



Jim S H [r) = - 



4noc 



r— >ro 



ro\ 



lim 6 ¥ {r) =0. (B3) 



This gives the standard Debye screening solution 

5u{r) 



—Area 
— =rrr exp 
\r - r \\ 



where the Debye screening length is 

dn[ji] 



A n = 4:71 0C- 



dp 



(B4) 



(B5) 



On distances larger than Ad, the charge of the electron 
is effectively screened. In particular, we can neglect the 
influence of the external electric field on motion of the 
electron if Ad is small compared to the de Broglie wave- 
length A = h/p of the electrons: 



1 « 



A 

Ad 



1 



mv 



I47T0C 



dn[ji] 



(B6) 



In the ultrarelativistic limit (A5), one has dn/dfi « 
}i 2 /tz 2 + T 2 /3 whereas in the Boltzmann limit (A10), 
one has dn/dji « Uq/T. 

The relevant electron velocity scale in our problem 
is v ~ 10 _3 c (T r>u 1 eV). This is the typical scale for 
electrons ionized from neutral hydrogen. For these elec- 
trons, one immediately sees that the screening becomes 
significant in the Boltzmann regime once the density is 
larger than 



m 2 v 2 T 

n > — ~ O.Iwb 

4tzoc 



(B7) 



= {mT/2n) 3/2 is the typical density in the Boltz- 
mann regime. The typical electric fields in this regime 
are E ~ V(p ~ eT/2z a , which yield an ionization po- 
tential of Ionize ~ eE^o ~ 0iTa^/7.z a rsj 1 meV <C 13 eV. 
Thus, the ionization of incoming hydrogen will not oc- 
cur until the density is much larger, by which point the 
screening will be highly efficient. Contrary to the argu- 
ments presented in [45], electrons depositive via neutral 
hydrogen can easily penetrate deeply into the electro- 
sphere producing the direct-annihilation emissions dis- 
cussed in Sec. V. 

Finally, we point out that, in the ultrarelativistic 
regime, we may express (B6) in terms of the Fermi mo- 
mentum: 



A 

Ad 



Pf 
V 



(B8) 



Thus, for highly energetic particles, screening is irrele- 
vant. The would apply to the ~ 5 MeV positrons ejected 
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from the core of the nuggets that we argued are respon- 
sible for the diffuse 10 keV emission [4]. These see the 
full electric field which is responsible for their deceler- 
ation and ultimately for the emission of their energy 
in the 10 keV band. This point is somewhat irrelevant, 



however, as the positron electrosphere of the antimat- 
ter nuggets cannot screen a positive charge (screening 
is only efficient for particles of the opposite charge than 
the ion constituents). 
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